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I . INTRODUCTION 

This report contains theoretical investigations on plasma processes 
in the Kaufman ion thruster. The temporal development of the neutraliza- 
tion of the ion beam by electron injection is analyzed in detail for dif- 
ferent physical conditions. The purpose of the analyses is i) to under- 
stand the formation of the steady-state, neutral ion beam, and ii) to 
predict theoretically the most favorable conditions for neutralization 
and the characteristic lengths and times required for neutralization. 

Among other things, the investigations indicate that the neutralization 
is brought about by nonlinear, discontinuous neutralization waves, and 
that a complete neutralization (without a periodic space charge structure) 
requires a strong intercomponent, collective momentum transfer between 
the collisionless electron and ion beams, i.e. and irreversible dissi- 
pation mechanism. Since the electron gas moves, relative to the ion 
beam, with supersonic drift velocity, a two-stream instability results 
which leads to a strong collective ("turbulent") momentum exchange be- 
tween the electron and ion components. Macroscopically , this type of 
intercomponent momentum transfer is described by Buneman's relaxation 
frequency in the electrohydrodynamic equations for the electron-ion beam 
system. This approach results in a quantitative, transient neutralization 
theory and a neutralization length which is compatible with experiment. 

CHAPTER II. The lateral neutralization of ion beams is treated by 
standard mathematical methods for first order, nonlinear partial differen- 
tial equations in the special case that the electrons are injected with a 
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density n Q equal to the ion density N. It is shown that the electron gas 
moves in form of a discontinuous step wave into the ion beam leaving the 
ion gas completely neutralized behind the head of the wave. The speed 
of the neutralization front is considerably increased by application of 
an auxiliary electric field in the direction opposite to the electron 
motion. 

CHAPTER III. By means of a von Mises transformation, a closed form 
analytical solution is derived for the transient, lateral beam neutrali- 
zation by electron injection with a density n Q different from the ion 
density N. It is shown that the electron gas moves in form of a discon- 
tinuous neutralization wave (1/2 < n /N < °°) or neutralization shock wave 

o 

(0 < n Q /N < 1/2) into the ion beam at the indicated density ratios. The 
neutralization wave produces a periodic over- and under-neutralization 
in the ion space traversed if ^ N at the injection plane. After its 
generation by the neutralization wave, the periodic space charge structure 
remains stationary if a dissipation mechanism is not considered. In this 
idealized, dissipation-free treatment of the transient neutralization 
problem, a complete neutralization is achieved only if the electrons are 
injected with a density n Q equal to the ion density N. 

CHAPTER IV. A nonlinear theory of the longitudinal ion beam neutral- 
ization is developed by means of the von Mises transformation in which inter- 
component momentum transfer by strong collective interaction of the Buneman 
type is considered as a dissipation mechanism. It is shown that the elec- 
tron gas moves in form of a discontinuous, but damped neutralization wave 
downstream the ion beam. As a result, the extent of the periodic over- and 
under-neutralization produced by the neutralization wave decreases within 
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a few relaxation lengths downstream of the neutralizing grid to zero. 

In a typical case, it is shown that a complete neutralization of the 
ion beam is obtained within a few centimeters downstream of the grid. 

CHAPTER V. By means of the Lenard-Balescu equation, the inter- 
component momentum transfer between stable, collisionless electron 
and ion components is calculated. The stable velocity distribution 
functions give rise to a momentum exchange between the electron and 
ion components which is too weak to explain the anomalous short neu- 
tralization lengths observed experimentally. For this reason, this 
kinetic theory is applied exclusively to the electrical current trans- 
port in stable, collisionless plasmas and a discussion of plasma 
conductivity. 

An investigation on neutralization shock waves has been initiated 
in the same research period but is not included in this report. This 
work will be communicated as soon as it has been completed. 

The investigations contained in this report represent preliminary 
results. The final (extended) versions of these investigations will be 
communicated through publications in the applied physics literature at 


a later date. 
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II. TRANSIENT ION NEUTRALIZATION BY ELECTRONS 


ABSTRACT 

The nonlinear initial-boundary-value problems describing the 
lateral neutralization of ion beams for the cases that 1) an 
auxiliary electric field accelerates the electrons into the ion 
space, and 2) the electrons are injected into the ions space at 
a prescribed current density are treated. Analytical solutions are 
derived which give the position and speed of the neutralization 
front in dependence of time, and the temporal development of the 
electron density, velocity and electric fields during the 


neutralization process. 
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The problem of the spreading and neutralization of rarified ion 

1 - 4 ) 

beams is encountered in electrostatic propulsion systems. The 

neutralization has been analyzed for the homogeneous ion beam model 

3 — 8 ) 

in the one-dimensional steady-state case. These investigations 

indicate that the electrons are distributed nonuniformly across the 

3 _ 8 ) 

beam and that the steady-state neutralization is incomplete. 

In the following, the nonlinear initial-boundary-value problem 
of the neutralization of a quasi-homogeneous ion beam is treated. 

The purpose of the investigation is to understand the most elementary 
physics of the temporal development of the neutralization process. 

Two simple models are considered for this purpose, which permit the 
study of the dependence of the neutralization on 1) the electron inject- 
ion mechanism and 2) the auxiliary electric field which moves the elect- 
rons of the neutralizer plasma into the ion space. 

3_8 

Following previous work on steady-state neutralization, a slab 

type ion beam of width a and infinite length is considered, which is 

3—8 

treated as homogeneous. The homogeneous ion beam model is not very 

realistic but required for mathematical reasons in order to render an 

analytical discussion feasable. The electrons are injected from the 

side into the beam, i.e. transverse to the beam velocity (Fig.l). This 

3-4 

arrangement is known as lateral neutralization . Only one-dimensional 
electron motions transverse to the beam velocity are analyzed, i.e. down- 
stream convection of electrons by the beam is disregarded. 
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THEORETICAL PRINCIPLES 


1 9) 

The electron gas moves in form of a nonlinear wave Ott < n /N < «)— 

l o 

or shock wave (0 £ h Q /N < y)— ' ^ into the ion gas at the indicated ratios 
of electron density YI q to ion density N . In order to achieve a 
complete neutralization (n = N) , the condition iX o = N has to be 

satisfied at the plane of electron injection. In this case, inertial 

effects play no role in the macroscopic electron motion, provided that 
the momentum relaxation time is small compared to the time it takes the 
electrons to penetrate the ion beam. 

With disregard of inertial effects, the field equations for the 

(* 

• 

drift velocity v and density, n of' the electrons and the self-" 
consistent electric field t are in presence of a homogeneous back- 
ground (N) of ions: 

V = (e/m) x E , (1) 




V*E = 4-rre(ft-N) , 


( 2 ) 

(3) 


when 


ID 


x = (M/m)"^3 2ir/u) , w = (4irNe^ /m)'*'^ 


(4) 


x is the collective relaxation time of the momentum density nmv and 
a) is the characteristic oscillation frequency of the electrons of 
elementary charge e < 0(m = electron mass, M = ion mass). 

Equation (1) is the simplified equation of electron motion which 
results in the absence of inertial effects for the collisionless 
electron-ion system [Eq. (4)]. A velocity field V is defined by 
Eq. (1) only in regions where E 4 0 , see Eqs. (2 ) — (3) . 
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A set of first order differential equations similar (but with a 

collisional relaxation time) has previously been used in the analysis 

12 ) 

of the collection of electrons in the pulsed ion chamber, — and 
electron conduction phenomena in solids The Eqs{l)-(4) are 
being applied herein to the transient, non-oscillatory (n Q = N) 
neutralization of ion beams by electrons. 
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NEUTRALIZATION IN AN ELECTRIC FIELD 

4 - 

As a model, an ion beam (V) of homogeneous density N is 
considered which is bounded in the planes x = 0 and X = a by grids 
which are impermeable except to electrons, Fig. 1. At the beginning 
of the neutralization, t = 0 , the adjacent space X < 0 is occupied 
by a homogeneous plasma reservoir of electron ion density, tt + = N , 
and a potential difference U = <f>(a) - 0 is quasi-instantaneously 

applied to the grids at x = 0 and x = a. The auxiliary electric field 
E o = -U/a represents a weak test field, i.e., the voltage U is assumed to 
be extremely small compared to the voltage equivalent of the beam energy: 

| e |u««— MV 2 = 10 3 - 10 4 e-volt . (5) 


Under the influence of the space charge field of the ions and 
the auxiliary field E q = -U/a , the electrons move through the grid 
plane X = 0 into the ion space, 0 a . In accordance with 

Eqs. (l)-(4), the resulting neutralization process is described by 
the nonlinear intial-boundary-value problem : 


9n 

at 



3E _ , 

— = a (n - V , 


where 


n(x, t = 0) = H(- x) , 
1 

/ E(x,t) dx = H(t) , 
o 

and 


( 6 ) 

(7) 

( 8 ) 
(9) 


a = 4teaN/(-U/a) > 0 


CIO) 
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In Eqs. (6)-(9), n(x, t) , E(x, t) , x and t are nondimensional 
variables which are related to the associated dimensional variables 
[Eqs. (1) — (4) ] by 

n = n/N , E = E/(-U/a) , (11) 

and 

x = x/a , t = t/ (-a/— - t) . (12) 

m a 

The reference constants in Eqs. (11)-(12) have an evident physical 
meaning, e.g., t Q = -a/(eUT/ma) is the transit time of an electron in 
the field -U/a for the distance a . 

Combining of Eqs. (6) and (7) yields for the derivative dn/dt 
along the stream line of an electron fluid element 

E -ante - 1) . (13) 

The solution of Eq. (13), which satisfies the initial condition in 
Eq. (8), is obviously: 

n(x, t) = H[s(t) - x] , (14) 

where 

s(t - 0) - 0 , (15) 

since 

[ ^ - E(x, t)] 6(s - x) = -aH(s - x)[H(s - x) - 1] 

for x = s(t) , (16) 

and 


E( s, t) 


(17) 
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by Eqs. (13)— (14 ) and Eq. (1), respectively. It represents s(t) the 
position coordinate of the neutralization front. It is seen from Eq. 

(14) that the initial condition, Eq. (8), determines the form n(x, t) 
of the solution for t >_ 0 . n(x, t) is a step function, which is 
discontinuous at x = s(t) , Fig. 2. 

Integration of Eq. (7) from x = s to x gives, under consideration 
of Eqs. (14) and (17) r 

E(x, t) = — + a(s - x)[ 1 -H(s - x)] , (18) 

at 

whence 

= H(t) + \ a(l - s) 2 (19) 

by substitution of Eq. (18) into the boundary condition in Eq. (9). The 
solution of Eq. (19), which satisfies the initial condition in Eq. (15), 
gives the position coordinate of the moving neutralization front: 

s(t) - 1 + (2/a) 1/2 tg[(a/2) 1/2 t -arctg(a/2) 1/2 ] ^0 , 

t 0 . 

(20) 

It follows for the speed ds/dt of the neutralization front in the 
external electric field: 

= cos 2 [ (a/2)‘*' //2 t - arctg(a/2)^ 2 ] >_ 1 , t >_ 0 . (21) 

The Eqs. (14) and (18), in which s(t) and ds/dt are given by Eqs. 

(20) and (21), constitute the nondimensional solutions for the fields 
n(x, t) , E (x, t) and v(x, t) . 

The time required for the neutralization front s(t) to penetrate 
the ion beam (0 x 1) is determined by s(t) = 1 . Accordingly, 
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the dimensionless neutralization time is by Eq. (20) 

t N = (2/a) 1/2 arc tg(a/2 ) 1/2 « ~ (2/a) 1/2 , (a/2) 1/2 » 1 . 

( 22 ) 

Since t^ = £^/£ Q » ^o ^ ^o ^ an< ^ a ^ ^ £ E 3 S * (10)-(14)], the 

dimensional neutralization time decreases with increasing 

E = -U/a . 
o 

In Fig. 3, the dimensionless position coordinate s(t) of the 
neutralization front is shown in dependence of the dimensionless time 
t with a as a parameter [Eq. (20)]. The line s(t) = 1 intersects 
with the s(t) curves at t = t^(a) , t^ being larger for smaller 
a-values. As to the dependence on the auxiliary field E q = -U/a , 
it is noted that a 'v* E q ^ and t £E q , i.e., that, for any given 
dimensional time £ the dimensional position coordinate 4 (£) is 
larger for larger E Q -values. 

In Fig. 4, the dimensionless velocity ds/dt of the neutralization 
front is shown in dependence of the dimensionless time t with a as 


a parameter. It is seen that, for any t < t N , ds/dt is larger for 

larger a-values. As to dependence on the auxiliary field E Q = -U/a , 

it is noted the dimensional velocity d4/d.£ = ( a/£ ) 'ds/dt increases 

o 

with increasing E q for any fixed instant £ = t£ Q , since £^ 'v E^ ^ 
and a ^ E Q . 

In Fig. 5, the dimensionless neutralization time t N is represented 

versus a [Eq. (22)]. The dimensional neutralization time = t^-t 

decreases with increasing E = -U/a , since £ E ^ and a 'u E^ ^ . 

o o o o 


As a numerical illustration, consider a mercury ion beam. In this case: 
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2 2 

-t = -] — r 77 - = 2.369 x 10 — N 

o |e| Ut U 

Hence 

a = 1.810 x 10 18 , x -1 = 1.249 x 10 10 sec -1 , t = 7.107 x 10 2 sec 

o 

16 —3 —6 

for a = 10 cm, N = 10 cm , and U = 10 /300 cgsu. Accordingly, 

the neutralization time is by Eq. (22): 

t N = 1.105 x 10" 9 , = 7.853 x 10 " 7 sec . 

It is seen that the neutralization time is rather short, although 
the driving test field is extremely weak. In this connection, it should 
be noted that n(x,t) = 1 for x _< s(t) [Eq. (14)], i.e., a net space 
charge field does not exist behind the neutralization front s (t) . 

Because of the disregard of inertial effects in the electron 
motion [Eq. (1)], the formulas derived in this section are only applicable 
to situations in which the neutralization time is large compared to 
the momentum relaxation time, >> x . 
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NEUTRALIZATION BY ELECTRON INJECTION 


As a model, an ion beam (V) of homogeneous density N is 
considered which is bounded in the planes x = 0 and X = a by grids 
which are impermeable except to electrons. Fig. 6. Through the plane 
X = 0 , electrons are injected quasi-ins tantaneously from the region 
X < 0 with a density = N , but an arbitrary current density 
/ = n Q eV Q at the beginning t = 0 of the neutralization (U = 0) . 

In accordance with Eqs. (l)-)4), the associated neutralization 
process is described by the . nonlinear initial -boundary -value problem: 


where 


8n 

at 


al (nE > 


3E 

3x 


mx(n-l) , 


(23) 

(24) 


n (x, t = 0) = H(- x) , (25) 

n(x = 0, t)v(x = 0, t) = H(t) , (26) 

and t , w is defined in Eq. (4). In Eqs. (23)-(26), n(x, t) , 
E(x, t) , x and t are nondimens ional variables which are related 
to the associated dimensional variables [Eqs. (1)— (4) ] by 

u 

n = n/N , v = u/v , E = E/<^-^) , (27) 

o e t 


and 


x = x/(u q / oj) , t = ut . (28) 

* 

The reference constants in Eqs. (27)-(28) have an eyident physical 

meaning, e.g. , E q = mv^/ez is the electric field that accelerates 

an electron to the speed V in x seconds. 

o 



Combining of Eqs. (23) and (24) yields for the derivative dn/dt 
along the streamline of an electron fluid element 

-^7 + E = -o)tn(n-l) . (29) 

The solution of Eq. (29), which satisfies the initial condition in Eq 
(25), is 


n(x, t) = H[cr(t) - x] , (30) 

where 

a(t = 0) = 0 , (31) 

since 

[^~ - E(x, t)3 d(a-x) = -arrH(a-x) [H(a-x) - 1] 
for x = a(t) , (32) 

and 

jjr=E(a, t) (33) 

by Eqs. (29)-(30) and Eq. (1), respectively. Equation (30) indicates 
that the initial condition, Eq . (25), determines the form of n(x, t) 
for t >_ 0 . The solution n(x, t) is a step function which is 
discontinuous at x = o(t) , the position of the neutralization front. 
Fig. 6. 

Integration of Eq. (24) from x = a to x gives, under 


consideration of Eqs. (30) and (33): 
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The position a(t) of the neutralization front is determined by the 
boundary condition in Eq. (26), which gives 

H(ct){^- + uto [l-H(a-x) } = H(t) (35) 

by Eqs . (30) and (34). It follows by integration of Eq. (35), since 
a (t) >_ 0 for t >_ 0 : 

a t 

/ H(a') da’ = / H(t’) dt’ . (36) 

0 0 

Hence, the position coordinate of the moving neutralization front is: 

a(t) = t , t 0 . (37) 

The speed da/dt of the neutralization front is invariant (constant 
electron injection velocity): 

% = 1 , t > 0 . (38) 

Cl t 

The Eqs. (30) and (34), in which a(t) and da/dt are given by Eqs. 
(37) and (38), constitute the nondimens ional solutions for the fields 
n(x, t) , E(x, t) and v(x, t). 

The time required for the neutralization front a(t) to penetrate 
the ion beam (0 x aw/v ) is determined by a(t) = aco/u Q [Eq. 
(28)]. Hence, the dimensionless neutralization time is by Eq. (37) 

C N = aaj/v o • (39) 

In Fig. 7, the (dimensionless) position coordinate a( t) and 
speed da/dt of the moving neutralization front are shown versus t 
in accordance with Eqs. (37) and (38), respectively. 
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The dimensional speed of the neutralization front and the dimensional 
neutralization time are by Eq. (38 ) — ( 39 ) and Eqs. (27 ) — (28 ) : 

V = j /Ne (40) 

o o 

and 

= a/V o (41) 

where V and j are the steady-state velocity and electrical current 
o o 

density with which the electrons are injected [Eq. (26)]. Accordingly, 

9 li" 0 l 7 

= 2.082 x 10 — — — = 6.246 x 10 cm/sec 

and 

t - 4.803 x 10~ 10 -S3- , 1.601 x 10" 7 sec 

hoi 

q n 2 O 7 

for a = 10 cm, N = 10 cm , and |/ | = 10 amp/cm = 3 x 10 cgsu. 

The speed of the neutralization front do/dt is identical with 

the speed u with which the electrons are injected (ft -N) which, in 
o o 

turn, determines the neutralization time In this case, the 

neutralization time can be shorter than the momentum relaxation time if 
the injection velocity is sufficiently large. 
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REMARKS 

Closed form analytical solutions have been obtained for the 
nonlinear initial-boundary-value problems describing the transverse 
neutralization of idealized ion beams in the cases that 1) a weak 
auxiliary electric field pulls the electrons into the ion space, and 
2) a constant electron current is injected through the neutralizer 
plane. In order to achieve a complete neutralization, the electron 
density at the neutralizer plane (x = 0) has been assumed to be equal 
to the ion density (N) . In the cases 1) and 2), the neutralization 
occurs through an electron flow with a discontinuous front, ahead of 
which the electron density is zero and behind of which the electron 
density is equal to the ion density (discontinuous solutions). 

In case 1), the speed of the neutralization front can be considerably 
enhanced by increasing the intensity of the auxiliary electric field. 

In case 2), the speed of the neutralization front is given by the 
velocity with which the electrons are injected (the electron gas 
behaves like an incompressible fluid). 

The main outcome of these simple theoretical considerations is that 

a short, relaxation-free neutralization is possible if the electrons are 

injected with a density equal to the ion density, = N. In the 

general case, n ; 4 N, the neutralization time is much longer because 

9 

of the occurrence of neutralization space charge waves , which exhibit 
only a weak damping. As a mathematical result, it is pointed out that 
the solutions of the nonlinear partial differential equations of first 
order are discontinuous. More realistic, continuous solutions are to 
be expected if effects associated with second order partial derivatives 
are taken into consideration, such as diffusion. 
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Electron Density Field n(x, t) Versus x at Time 


a»IO a = 10 a=!09 8 7 



FIG. 3: Nondimens ional Position Coordinate s(t) of the Neutralization Front 








FIG. 5: Nondimensional Neutralization Time t„ Versus Parameter 
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FIG. 6: Discontinuous Electron Density Field n(x, t) Versus x at Time 
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III. NONLINEAR ELECTRON-ION NEUTRALIZATION WAVES WITHOUT 
DISSIPATION (LATERAL INJECTION) 

ABSTRACT 

The nonlinear boundary-value problem describing the dynamics of 

an electron gas, which is injected with a prescribed current density 

from a plane (emitter) into an ion gas, is solved in closed form by 

means of a von Mises transformation. The electrons and ions are 

assumed to interact predominantly through the self consistent space 

charge field, while non-conservative forces are disregarded. It is 

shown that the electron gas propagates in form of a regular nonlinear 

wave < n /N < °°) or shock wave (0 < n /N < 4) into the ion gas 
o o £ 

at the indicated ratios of electron injection (n ) to ion (N) 

o 

density. Graphs of the electron density, velocity, and electric space 
charge field in the neutralization wave are presented for various times, 
as well as plots of the time-dependence of the position coordinate 
of the neutralization front. The neutralization waves lead to a com- 
plete neutralization i) if n^ = N at the injection plane, and ii) to 

a spatially periodic over- and under-neutralization if n / N . 

o 

These results are strictly valid only for times which are small com- 
pared to the relaxation time for the dissipative momentum transfer 
between electrons and ions . 
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The formation of an electrically neutral plasma by injection of 
electrons into a rarefield ion gas is a nonlinear process in which 
electron-ion interactions through the self-consistent space charge 
field play a predominant role. In the final satages of the neutraliza- 
tion, diffusion and dissipative, intercomponent momentum transfer become 
significant as dispersion mechanism for periodic, local neutralization 
unbalances . 

As will be shown, the electron gas moves into a collisionless ion 
gas in form of a nonlinear wave (% < n Q /N < «•) or shock wave 
(0 < n Q /N < %) at the indicated ratios of electron (n Q ) and ion (N) 
densities. For this reason, it is distinguished between i) "nonlinear 
neutralization waves" (no overtaking of particles) and ii) "neutraliza- 
tion shock waves." Both phenomena exhibit a discontinuous neutralization 
front (ahead of which the electron density is zero) , which moves with 
the "neutralization speed" into the ion space. 

The nonlinear theory to be presented is concerned exclusively 
with ordinary neutralization waves (^ < n Q /N < <*> at the injection 
plane). The ions are assumed to be distributed homogeneously at the 
beginning of the neutralization (initial condition) . The redistribution 
of the ions (Coulomb repulsion) during the period which the electrons 
need to penetrate a distance of the order of several neutralization 
wave lengths (x q ) in the ion space is commonly negligible because 
of the large inertia of the ions (nu >> m e ). 
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Further, as the electrons proceed in neutralizing the ion gas, 

the Coulomb repulsion is more and more weakened. It is recognized 

that the ion gas can be treated as a homogeneous space charge 

background for times which are large compared to the characteristic 

period (t^) of the neutralization process. 

The concept of a homogeneous ion gas with uniform drift velocity 

bounded by a plane plate has been proposed as a model for an ion 
1-2 

beam. In this respect, the neutralization wave theory is 

applicable to the transient neutralization of ion beams of electro- 

3-4 

static propulsion systems. The idealized theory should provide 
a qualitative understanding of the neutralization of ion beams in 
free space. Consideration of the inhomogeneity and spreading of free 
ion beams does not seem to be mathematically feasible except by means 
of a considerable numerical expenditure. The previous work has 
been concerned with the steady-state neutralization of quasi- 
homogeneous ion beams . ^ ^ ^ ^ 
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THEORETICAL PRINCIPLES 


Subject of the considerations is a collisionless electron gas 
which is penetrating a quasi-homogeneous ion gas of density N . The 
ion gas is limited to the slab space 0 .<_ x < °°(|y| < 00 , |z| < 00 ) . 

The plane at x=0 serves as an emitter which injects electrons at a 

! .t ' V - ‘ t 

prescribed current density j(t) as indicated in Fig. 1. The field 
equations describing the density n(x,t) , velocity v(x,t) and 
longitudinal space charge field E(x,t) are^ (electron charge = 
e < 0 , electron mass = m) : 

3v " e 


3v 

3t 


tt - + v — = — E , 
3x m 


( 1 ) 


3n , 3 , x n 

a? + 57 (nv) ” 0 • 


^ - -4ire(N - n) 

OX 


( 2 ) 

(3) 


The injection of the electrons through the plate at x=0 with the 
current density j(t) = en(0,t)v(0, t) is taken into consideration 

i f 

by means of the boundary conditions: 


[n(x,t)v(x ,t)] x=Q = j (t) / e > 0 


(4) 


n(X,t) K=0 * ”o 


(5) 


The continuity of the displacement and convection currents across 
the interface at x=0 requires 

^St^ + en(x,t)v(x,t)] x=0 = j * (6) 

The Eqs. (4)- (6) imply for the displacement current the boundary 
condition, 3E(x,t)/3t = 0 for x = +0 . 
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\ 


In Eq. (1), the pressure gradient has been neglected in comparison 
to the electric field force (analogous to the zero-temperature 
plasma approximation). The Eqs. (1) — (2 ) can be linearized by either 

g 

transforming to Lagrangian coordinates or by introducing a stream 

9 

function in accordance with the von Mises transformation. Both 
methods have been applied in the nonlinear theory of Langmuir 
oscillations of plasmas of infinite extension^ ^ and can be 
extended to the analysis of nonlinear phenomena in non-neutral 
particle gases of finite extension with proper boundary conditions. 

The nonlinear boundary-value problem defined in Eqs. 

(l)-(6) is subject to the von Mises transformation. Let a stream 

function ^ = ^(x,t) be introduced, which automatically satisfies 

/ 

Eq. (2), by the relations 


dip _ nv dip _ n 
3t ° N ’ 3x N ' 


(7) 


The stream function is constant along the trajectories of the electron 
fluid since 


d^ _ a^ dx jhp 

dt 3x dt at 


x 


n 

- N V 


nv 

N 


= 0 


( 8 ) 


If the Jacobian is J{ (x,t)/(\fj,t) } = [3x/a^] t ^ 0 , it is possible 
to change from the old (x,t) to the new , t ) independent 

variables : 

x = x(^,t) , t = t ; 

n(x,t) -*■ n(i|>,t) , E(x,t) -*■ E(if>,t) , v(x,t) -* v(4>,t) . (9) 


Since 


dx __ 

ax 

dip , ax 


_ M 

” dt 

3ip 

t dt at 

ip 

at 


v 
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3v 

at 


<P 


3 2 


x 


3t 


and 


dv _ 3v 

3v 

+ v 

3v 

dip 3v 

3v 

dt = 3t 

3x 

3i|t 

dt 3t 

, 3t 


X 

t 

t 

♦ 


it follows that 


3v . 3v 

3l + 


3 2 x 

2 

3t 


( 10 ) 


(ID 


( 12 ) 


Thus, by changing to the independent variables 0Jt,t) and 
considering the interrelations in Eqs. (8), (9) and (12), the nonlinear 
Eqs. (1) — (3) are reduced to the linear system: 



e 3E 
m 3x 


- <o 2 (l +ff) 


(13) 

(14) 


where 

to = (4TrNe 2 /m) 2 (15) 

is the electron plasma frequency. Integration of Eq. (14) from x=0 
to x gives 


^ E(ijt,t) = —to 2 + x - [ip (0 , t ) + E(0,t) ] } 

mto 

Integration of Eqs. (2) and (3) from x=0 to x gives, under 
consideration of Eq. (7), an equation physically equivalent to 


Eq. (14), 


i- M + eN 11 

4 it 3t 3t 


j(t) , 


(16) 


(17) 


where the integration "constant" is j (t) by Eq. (6). Integration 
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of Eq. (17) from t=0 to t yields an alternate form of Eq. (16), 

^ E(<fj,t) = -w 2 {<|i - O(x,0) + -^2 E(x,0)] - — / j(t')dt') . 

mo) 0 

( 18 ) 

A comparison of Eqs. (16) and (18) reveals that the various integration 
"constants," j|j(x,0) , . . . ,E(0,t) , are related in a simple way to 
x and j(t): 

<Kx,0) + E(x,0) = -x , (19) 

mco 

e i , c 

♦ (0,t) + — ^ E(0,t) = ^ / j(t')dt' . (20) 

mo) 0 

Elimination of E (ip , t ) from Eq. (13) by means of Eq. (16) 
or Eq. (18) reduces the nonlinear system in Eqs. (l)-(3) to an 
inhomogeneous, linear oscillator equation for x = x(i|j,t): 

+ od 2 x = -w 2 i|j + oj 2 f (t) , (21) 

9t 

where 

t 

f(t) = I j(t')dt* . , (22) 

0 


The general solution of Eq. (21) is by Lagrange's method of the form 


where 


x(i|),t) = A(^)sin wt + B(\|0cos wt - ip + F(t) 


t 

F(t) = oj / f(t')sin a)(t-t')dt' 
0 


(23) 


(24) 


with F(0) = 0 and F'(0) = 0 . The functions A(iJj) and B(ijO are 
integration "constants" which are determined by the respective 
initial and/or boundary conditions for the fields n(x,t) and 
v(x,t) , which are related to ijj(0,t) and \p(x,0) by Eq. (7): 
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t 

\p(0,t) - <p(0,0) = - / [n(x , t ' )v (x, t ' ) ] _ dt' , (25) 

0 

and 

1 x 

iJj(x, 0) - ^ (0 , 0) = — / [n(x* ,t)] t=0 dx' . (26) 

By means of the general solution for x(ij>,t) in Eq. (23), one 
obtains E(>|j,t) by substitution of x(i^,t) into Eq. (16), v(i|i,t) 
by partial differentiation 3/3t)^ of x(^,t), and n(ifj,t) by 
partial differentiation 3/3<|O t of x (i^ , t ) : 

E(^,t) = - — w 2 [A(i|/)sin cot + B(i|j)cos ait - f(t) + F(t)] , 

6 (27) 

v(i(j,t) = oi[A(^)cos tot - B(i|»)sin cut] + F'(t) , (28) 

n(\p,t) = N/[l - dA ^|^ sin wt - cos ut] . (29) 


The Eqs. (27)— (29) represent a parametric solution for the fields 
E(x,t) , v(x,t) and n(x,t) in terms of >p and t . The function 
\l> = i|/(x,t) is given implicitly by Eq. (23). 
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NEUTRALIZATION WAVES 

The boundary-value problem describing the motion of the electron 

gas in the ion space x ^ 0 is defined in Eqs. (1) — (6 ) . In order to 

keep the mathematical expenditure a minimum, an electron current injection 
corresponding to a step-impulse is considered: 

j(t) = i Q H(t) , (30) 

where 

H(t) =1, t > +0 , (31) 

and = 0, t £ -0 , 

and 

i = en v <0 , (32) 

o o o 

so that 

a = Ne/ i > 0 . (33) 

Eq. (30) represents a quasi-instantaneous injection of the current 
density i Q . 

Since n(0,t)V(0,t) = j(t)/e is a known boundary-value [Eq. (30)], 
Eq. (25) is relevant which gives 

<K0,t) - <K0,0) = o" 1 H(t ) t = *(0,0 +— E(0,t) (34) 

2 

mu)^ 


by Eq. (20). Comparison of these relations indicates that E(0,t) 
is constant, 


— E(0,t) = *(0,0) 5 * , (35) 

mw 2 0 

and 3E/3t =0 at x=0 . In the absence of surface charges at 
x = +0 , it is = 0 and E(0,t) = 0 . Equation (34) is a relation 
between ip and t at x = 0 , 

ip - ip = a 1 H(t) t ; t = t , $ = ip (0 , t ) . (36) 

Evaluation of Eqs . (22) and (24) for the current density in 


Eq. (30) yields 

f(t) = a -1 H(t)t , (37) 

F(t) = — H(t) (wt - sin wt) . (38) 

ojo 

The boundary conditions in Eqs. (4)-(5) give v(x=0,t) = v(ijj,t) , 
and upon application to Eq. (28), 

A(i|j)cos wt - B(ijj)sin a)t = — [v(x=0,t) - F'(t)] , (39) 

W 

while 

A(tji)sin ut + B((ji)cos wt = * - F(t) (40) 


by Eq. (23) for x = 0 . One obtains by elimination from Eqs. (39)-(40), 
under consideration of v(x=0,t) = (N/n o )H(t)/o and Eq. (38): 

A (ip ) = ip sin wt + — [1 - (1- — )cos wx]H(t) , (41) 

O UiO ti 

o 

B (ip) = ip COS WT + — [(1- — )sin wt]H(t) , (42) 

o wo n 


where 
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since the structure of a function does not depend on the designation 
of the independent variable (tjj ■ ip , t ->- t) . 

Upon substitution for A(^) , B(^) , and f(t) , F(t) , in 

accordance with Eqs. (37 ) — ( 38 ) and Eqs. (41)-(42), respectively, one 
obtains from Eqs. (27 ) — (29) the parametric solutions: 


- = -(1- — )H(x)sin to (t — t ) + [H(x)-H(t) ]sin cot + wa\p cos co(t-x) , 
-mco / ea n o 

° (44) 


= H(t) -(1- — )H(x)cos co(t-x) + [H(x)-H(t) ]cos cot - wo\p sin co(t-x) , 
l/o n o 

o 

(45) 


nOJst) _ H(x ) ^ 

N H(x)-(1- — )H(x)cos co (t-x ) - — " —- -s in cot-coo^ sin co(t-x) 
n n co o 

o o 


(46) 


where 

N 

coax(^,t) = -uicnp + H(t)cot - (1 )H(x)sin co(t-x) 

n 

o 

+ [H(x) - H(t)]sin cot + cocij^cos co(t-x) , (47) 

and x = x (^ ) is given in Eq. (43). Equation (47) gives ip = ip(x,t) 
for every point (x,t) , where t = x for x = 0 [Eq. (35)]. 

For the purpose of discussion, it is suitable to introduce non- 
dimensional variables and fields in accordance with the substitutions 


x/x x , t/t -* t , ip/x -*■ \p ; 

° ° ° (48) 

E/E -*■ E , v/V ->■ v , n/N -*• n ; 

o ’ o * o 

x = 1/coo , t = 1/co , E = -mco/eo , V = l/o , N = N 
o o o o o 


If in addition surface charges are absent (i[> = 0), the Eqs. (44) — (47) 

become, for t >_ 0 and x >_ 0 : 
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E(ty ,t) = -(l-'^-)sin(t-ip) , 

o 

v(ip,t) = 1 -(1- ~ ) cos (t— if) ) , 

o 

n(ij),t) = 1/ [ 1 -(1- ~ ) cos ( t — <5 (if)) sin t] , 


(49) 

(50) 
(5L) 


where 

x(i|),t) = (t-i|0 - (1- ^-)sin(t-if)) , ip > 0 , (52) 

o 

and t has been eliminated by means of Eq. (43) for t > 0 , i.e., 
t ^ 0 and if) ^ 0 (4> q = 0) . Note that in Eqs. (49) -(52), N and n Q 
are the original, dimensionless densities. 

The Eqs. (49) -(52) describe a regular nonlinear wave if 
J = 9x(i|),t)/3if) 4 0 [n(if),t) ^ “ in Eq. (51)], i.e., if 

0 < N/n < 2 , or: h < n /N < 00 . (53) 

o o 

The position x = x(t) of the head of the wave (neutralization front) 

/V 

obtains from Eq. (47) as the limit x = x(ip -*■ +0,t) , 

x (t ) = t - (1- f-)sin t > 0 , t > 0 . (54) 

o 

Accordingly, x = kir for t = kn, k = 0 , 1, 2, 3,..., while 
x > t depending on t =/ ku and N/n Q > 1 . The speed of the 
neutralization front is 


dx(t) _ 
dt 


n N ^ 

-(1- — )cos 

o 


t = v(lj/,t) 


i |)=0 ‘ 


(55) 


It is seen that x(t) 0 and dx(t)/dt > 0 for t >_ 0 since 

0 < N/n <2 by Eq. (53). 
o 


In Figs. 2, 3, and 4 the (nondimens ional) fields n(x,t) , 

E(x,t) and v(x,t) are shown versus x with 0 £ t 10tt and 

= N/n = 0.1 and 1.1 as parameters based on Eqs. (49)-(51). At 
o 


o 
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a given time t = t , the neutralization wave has penetrated the 

A A A 

region 0 <_ x <_ x(t) , where x(t) is the position of the neutralization 
front at that time. The broken curve sections in Figs. 2-4 represent, 
e.g., the neutralization front at time t = 2tt , where x(2ir) = 2ir 
by Eq. (54). In Fig. 5, the position x(t) of the neutralization 
front is shown for arbitrary times t with a as a parameter 
[Eq. (54)]. 

In the region ahead of the neutralization front, x > x(t) , it is 
x < 0 and, hence, ip = 0 (\p^ = 0) or \p = ip^ (ip ^ 0) by Eq. (43). 
Accordingly, the fields in the unneutralized ion space are by Eqs. (44-46) 

n(x, t) = v(x,t) = 0 , E(x, t) = x-t , x > x(t) . (56) 

The Eqs. (49)-(52) demonstrate (see Figs. 2-4) that the 
neutralization wave produces an incomplete neutralization of the ions 
in the region 0 < x < x(t) , if a = N/n Q 4 1 . In periodically 
adjacent regions, an over- and under-neutralization is achieved (Fig. 2). 
The association space charge field changes periodically its direction 
(Fig. 3), while the velocity fields fluctuates periodically between 
the values v + = 1 ± [1 -(N/n Q )] , Fig. 4. If a = 1 , the 
neutralization wave reduces to a step wave with a neutralization speed 
dx(t)/dt = 1 [Eqs. (49)-(62) with N/n^ =1]. In the latter case, 
a complete neutralization is observed [n(x,t) = 1 for 0 j^x <^x(t))]. 

The reference values x , t , E , V , and N defined 

o o o o o 

in Eq. (48) determine the actual magnitude of the variables. As a 
numerical example, consider an ion gas (N) and current (i Q ) 
density given by 

9-3 - -2 -2 

N = 10 cm , i = 10 amp cm 
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Hence : 

x = 3.691 x 10 4 li |n" 3/2 = 3.502 x 10" 2 cm 
o 'o' 

t = 1.773 x 10 -5 N _1/2 = 5.607 x 10~ 10 sec 
o 

and 

E o = 2.227 x 10" 4 j i Q [N~ 1/2 = 2.113 x 10 _1 cgsu = 6.339 x 10 1 

volt cm 1 , 

V = 2.082 x 10 9 li IN -1 = 6.246 x 10 7 cm sec" 1 , 

O 1 0 1 

9 -3 

N = N = 10 cm . 
o 

The values of x q and t give the magnitude of the space and time 
periods of the neutralization wave (V q = x o /t Q ) . The (positive) 
amplitudes E , v + , and n + of the dimensional fields are by 
Eqs . (49)-(51) : 

E= | (1- ~ -)E o | , v ± - | [1 ± (l-f-)]V o | , n + = |N/t ± d-f-)]| 
o o o 

(57). 

In this idealized treatment of neutralization waves, dissipation 
mechanisms have been disregarded. Consideration of dissipation would lead 
to a damping of the wave amplitudes and thus to a dispersion of the local 
space charge unbalances, i.e. to a complete neutralization. The most impor- 
tant dissipation is caused by irreversible momentum exchange between the 
electron and ion components. By means of a dimensional argument, one finds 
from the inter component momentum relaxation time x and the characteristic 
parameters i Q , N, e as order-of-magnitude of the relaxation length for 
neutralization 

li /Ne 
o 


A » 


T 


(58) 
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FIG. 2: Electron density field n(x,t) versus x at time t = 10u for a = N/n 
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FIG. 3: Electric field E(x,t) versus 0 < x < x(t) at time t = IOtt for a = N/n 
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FIG. 4: Electron velocity field v(x,t) versus x at time t = 107T for a = N/n 
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FIG. 5: Position coordinate x(t) of the neutralization 
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IV. NONLINEAR THEORY QFt ELECTRON NEUTRALIZATION WAVES IN 

ION BEAMS WITH DISSIPATION 
(LONGITUDINAL INJECTION) 


ABSTRACT 

An analytical theory of nonlinear neutralization waves generated 
by injection of electrons from a grid in direction of a homogeneous 
ion beam of uniform velocity and infinite extension is presented. 

The electrons are assumed to interact with the ions through the self- 
consistent space charge field and by strong collective interactions, 
while diffusion in the pressure gradient is disregarded (zero- 
temperature approximation) . The associated nonlinear boundary-value 
problem is solved in closed form by means of a von Mises transforma- 
tion. It is shown that the electron gas moves into the ion space in 
form of a discontinuous neutralization wave, which exhibits a periodic 
field structure (incomplete neutralization) . This periodic wave 
structure is damped out by intercomponent momentum transfer, i.e., 
after a few relaxation lengths a quasi-neutral plasma results. The 
relaxation scale in space agrees with neutralization experiments of 
rarefied ion beams, if the collective momentum transfer between the 
electron and ion streams is assumed to be of the Buneman type. 
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A beam of positive ions attracts electrons from the surroundings , 

such as system walls and rest gases of an incomplete vacuum. Thus, a 

partial neutralization by electrons is a quite common phenomenon in 

1-2 

ion beams. In many technical applications, one is interested in 

producing a complete neutralization of ion beams by electron injection 

in order to reduce beam spreading and beam current limitations by 

3-4 

build-up of large space charge potentials. 

The neutralization of a rarified ion beam by injection of an 
electron-gas is a nonlinear process in which electron-ion interactions 
through the self-consistent space charge field play a predominant role. 
In the final stages of the neutralization intercomponent momentum trans- 
fer and to some extent also diffusion resulting from pressure gradients 
become significant as. despersion mechanisms of periodic, local neutral- 
ization unbalances. As will be shown, the electron gas moves into a 

collisionless ion beam in form of a nonlinear wave (1/2 5 n /N < °°) 

o 

or shock wave (0 < n /N < 1/2) at the indicated ratios of electron 

o 

injection (n Q ) and i°n (N) densities. For this reason, it is 
distinguished between i) "nonlinear neutralization waves" (no over- 
taking of particles) and ii) "neutralization shock waves." Both . 
phenomena exhibit a discontinuous neutralization front (ahead of which 
the electron density is zero), which moves with the "neutralization 
speed" into the ion space. 

As a model, a quasi-homogeneous ion beam if infinite extension is 
considered, i.e. the density N and velocity V of the beam are 
assumed to be uniform. The electrons are injected with a prescribed 
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current density j(t) in the downstream direction (j | |V) . The 
homogeneous ion beam model may be justified as i) a steady state 
condition, and ii) an initial condition for the following reasons: 

i) Since the beam is of infinite extension, a steady state exists 


with a homogeneous ion density N and a uniform flow velocity V . 
Because of the relatively large inertia of the ions (M >> m) , the 

perturbation of the homogeneity of the ion gas by injection of an 
electron stream from a grid is negligible for a period which is large 
compared to the characteristic time (t Q ) of the neutralization wave. 

ii) The redistribution of the ions during the period which the 
electrons need to penetrate a distance of the order of several neutral- 
ization wave lengths (x^) in the ion space is negligible because of the 
inertia of the ions (M >> m) . Further, as the electrons proceed in 
neutralizing the ion gas, the perturbing interactions through the 
self-consistent field are reduced. 

The homogeneous ion beam model is adequate for analyzing and 
understanding the basic properties of transient neutralization waves. 

The neutralization in a finite, inhomogeneous ion beam does not seem to 
be treatable by analytical methods . The previous analytical work has 
been concerned with steady-state electron-ion neutralization without 
dissipation in homogeneous ion beams, 
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NONLINEAR - BOUNDARY - VALUE PROBLEM 

An ion gas of homogeneous density N and uniform velocity V 
parallel to the x - axis is considered in the space |x| <_°° , |y| <_«> , 
|z| 00 . This infinite ion beam is bisected by a permeable grid in 

the plane x = 0 , from which electrons are ejected with the current 
density j(t) in the direction of V (Fig. 1). The field equations 
for the electron velocity v(x,t) , electron density n(x,t) , and the 
selfconsistent electric field E(x,t) are^: 


9v 

9t 


+ v 


3v 

9x 


- - E - v(v-V) , 
m 


(i) 


9n 

at 



9E 

9x 


4ire(N-n) 


( 2 ) 

(3) 


where e > 0 is the elementary charge and m the electron mass . 

The term -v(v-V) , describes the momentum transfer between the elec- 
trons and ions (v = relaxation frequency). The injection of electrons 
with the current density j(t) into the ion beam from the grid plane 
at x = 0 is taken into consideration by the boundary conditions: 


[n(x,t) v(x,t)] x=Q = j(t)/(-e) , (4) 

n(x,t) x=0 = % ’ (5) 

The semispace x < 0 is assumed to be shielded from electric dis- 
placement currents by the grid at x = 0 . Accordingly, the continuity 
of displacement and convection currents across the plane x = 0 gives 

''<*■« l x -0 * 


3 (t) . 


( 6 ) 
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The Eqs . (4) - (6) imply that the displacement current satisfies the 
boundary condition , 3E(x,t)/3t = 0 for x = + 0 . 

The electron and ion gases are assumed to be rarefied 
(n, N < 10 32 cm 3 ) and of low temperature (T < 10^°K) . The effects 
of binary collisions and linear electron-wave interactions are not 
discussed. Possible strong collective ("turbulence") momentum exchange 
between the electron and ion streams is discussed based on Buneman's 
relaxation frequency 3 

v = (m/M) 1/3 u)/2TT , (7) 

where 

ai = (47rNe 2 /m) 1/2 . (8) 

The Eqs. (1) - (6) describe the nonlinear-boundary-value problem 
of the transient neutralization process in the ion beam which takes 
place in the region x >_ 0 . Thermal effects are not taken into 


consideration . 
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METHOD OF SOLUTION 

The nonlinear-boundary-value problem defined in Eqs . (1) - , 

9-10 

(6) is solvable in closed form by a von Mises transformation 

Let a stream function t | i = i^(x,t) be introduced, which automatically 

satisfies Eq. (2), by the relations 


3\|) _ nv 3 jJj_ n 

3t ~ N ’ 3x N * 



The stream function is constant along the trajectories of the electron 
fluid since 


dip _ dx dip_ 

dt 3x dt 3t 

The Jacobian is J{ [ (x,t) / (i)j ,t) ] } = 3x/3i|>| t f 0 
For such solutions, it is permitted to introduce 
variables in place of (x,t) : 


n , nv 

= “F + r 


= o 


( 10 ) 

for n(x,t) < 00 . 

(ip,t) as independent 


x = x(ip,t) , t = t , 


n(x,t) + n(t|j,t) , E(x , t) -»■ E (ip , t) , v(x,t) -+v(ip,t). (11) 


Since 


v = 


dx 

dt 


3x 

3i}) 


dip , 3x 

dt + 


dx 

3t 


i.e. 


3v 

3t 


= jTjc 

ip dt 2 


( 12 ) 


and 


dv 

dt 


3v 

3t 


8v 

+ v — — 


3v 

dif; 3v 


dv 

dx 

X 

t 

3i|< 

dt at 

<P 

dt 


( 13 ) 



-52- 


it follows that 


3v 3v _ 3^x 

at v 3x = at 2 • 


( 14 ) 


Thus, by changing to the independent variables (ij>,t) and 
considering the interrelations in Eqs . (10), (11), and (14), the 
nonlinear Eqs. (1) — (3) are reduced to the linear system: 



S E ■ v( !t ■ v) 


(15) 


e 3E 
m 8x 


to 2 (1 + 


8K 

ax'* 


Integration of Eq. (16) from x = 0 to x gives 


(16) 


- E(*,t) - co 2 {* + x - [iK0,t) - “S’ E(0,t)]) . (17) 

m ^ 

mw 

Integration of Eqs. (2) and (3) from x = 0 to x gives, under 
consideration of Eq. (6), an equation physically equivalent to Eq. 
(16), 


1 9E 

4ir 3t 


eN lf 


j(t) , 


(18) 


where the integration "constant" is j(t) by Eq. (4). 
Eq. (18) from t = 0 to t yields another form of Eq. 

9 , t 

— E(^ ,t) = u) {ip - [^(x,0) 2 E (x,0)] + ^ / 

m mo) 0 


Integration of 
(17), 

j(t')dt’} . 

(19) 


A comparison of Eqs. (17) and (19) reveals that the various integration 
"constants," t|)(x,0), ... , E(0,t) , are related in a simple way to 


x and j(t) : 
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i|/(x,0) 2 E ( x >°) = “ x , 


mu) 


iKO,t) - -^2 E(0,t) = /j(t')dt' 


( 20 ) 


( 21 ) 


mu) 


Elimination of E(tp,t) from Eq. (15) by means of Eq. (17) or 
Eq. (19) reduces the nonlinear system in Eqs . (l)-(3) to an inhomo- 
geneous, damped, linear oscillator equation for x = x(ij),t) : 


3 2 x , 3x , 2 

— 7 ; + v t— + u) x 

3t 3t 


= vV + u) [f(t) - <J)] 


( 22 ) 


where 


f<t> = - s; / . 


(23) 


By Lagrange's method, the general solution of Eq. (23) is for 
u) > v/2 [Eq. (7)]: 

x(<J),t) = e [A(ij))sinfit + B (ip ) cosfit] - tj) + F(t) (24) 


where 


n = [u) 2 - (v/2) 2 ]^ > 0 , 


(25) 


and F(t) is a particular integral of the inhomogeneous differential 
equation. 


d 2 F dF 2 

— 9 + V -T- + U) F 

dt 2 dt 


= vV + u) f(t) , 


(26) 


i.e. 



-54- 


F(t) = - £- / f(t , )e v(t ’" t)/2 sinfi(t , -t)dt' 

u Q 

4- ^-[1 - e sinOt + cosftt)] , (27) 

w 

with F(0) = 0 and F'(0) = 0 . A(ifi) and B(^) are arbitrary 
functions which are determined by the respective initial and/or 
boundary conditions for n(x,t) and v(x,t) . By Eq. (9) 

■, t 

i|j(0,t) - <K0,0) = tt I [n(x,t')v(x,t’)] dt' , (28) 

0 

and 


ij;(x,0) - i(j( 0,0) = -| / [n(x* ,t)] t=0 dx' . (29) 

interrelates the boundary (x=0) and initial (t=0) values of \p(x,t) 
to those of n(x,t) and v(x,t) . 

By means of the general solution for x(i|j,t) in Eq . (24), one 
obtains E (ip , t ) by substitution of x(ip,t) into Eq. (17), v(^,t) 
by partial differentiation d/dt)^ of x(tj>,t) [Eq. (12)], and n(ip,t) 
by partial differentiation d/dtp)^ of x(ip,t) [Eq. (9)] : 

2 

E (i[» , t) = {e Vt / 2 [A(i|/)sinftt 4- B(ijj)cosfit] - f(t) + F(t)} , 

(30) 

v(^,t) = e [A(^) cosfit - B(40sinfit] 

- ^-[A(t}))sinnt + B(i^)cosfit]} + ^ , (31) 

l at 

n(i|j,t) = N{ 1 - e Vt ^ 2 [ sinfit + cosftt]} ^ . (32) 

dip 
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The Eqs. (30) - (32) represent a parametric solution for the fields 

,rv_ 

E(x,t) , v(x,t) , and n(x,t) in terms of ip and t . The function 
i|> = i}j(x,t) is given implicitly by Eq. (24). 
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LONGITUDINAL NEUTRALIZATION WAVES 

The theoretical investigation of neutralization waves generated 
by electron injection into a homogeneous ion beam requires specification 
of the time dependence of j(t). For mathematical convenience, an 
electron injection from the grid plane x = 0 in the form of step- 
impulse is considered: 

j(t) = i o H(t) , t S 0 , (33) 

where 

H(t) = 1, t _> +0 , 

= 0, t < -0 , (34) 

and 

i = - en v <0 . (35) 

o o o 

Hence , 

a = — Ne/i > 0 (36) 

o 

This model corresponds to a quasi-instantaneous electron injection 
with the current density i Q . 

A comparison of the Eqs. (21) and (28) indicates that E(0,t) is 
a constant independent of j(t), 

(e/mw 2 )E(0,t) = ip (0 , 0 ) = ip^ 

In absence of surface charges at x = 0 , it is E(0,t) = 0 and 
ip = 0 . 

1 A 
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The boundary-value n(x=0,t)v(x=0,t) = -j(t)/e is known through 
Eqs. (33) and (35). Accordingly, Eq. (28) is relevant which gives 

<Kt) - if* 0 = o 1 f(t), ip(t) = ^(x=0,t) . (37) 

The boundary conditions in Eqs. (4) - (5) give v(x-0, t)= -j(t)/n Q e , 
and upon application to Eq. (31) 

A(iJi)cosftt - B(iJj)sinftt = fi ■'■ e Vt ^ 2 { - F(t)] - — F'(t)} 

Z n e 

o 

(38) 

since 

A(\|))sinftt + B(i[))cos^t = e Vt/ ^ 2 [\}> - F(t)] (39) 

by Eq. (24) for x = 0 . The Eqs. (38) - (39) represent two independent 
relations from which one obtains by elimination: 

A ( 41 ) = e VT ^ 2 l> - F(t) ]sinfix 

+ fi _1 e VT/2 { y[i|j - F(t) ] - - F' (x) }cosfix , (40) 

o 

and 

B(tJj) = e VT ^ 2 [t|; - F(t)]cos fix 

-^ _1 e VT/2 { fa - F(x) ] - - F' (x) }sinfix (41) 

z n e 

o 

where 

o(\p - ^ q ) = f(x) , (42) 

since the structure of a function does not depend on the designation of 
the independent variable (i|» -*■ ip) . In Eqs . (40) - (41) , x is a 
function of ip defined by Eq. (42). 
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With the current density j(t) given by Eq. (33), the associated 
functions f(t) and F(t) in Eqs. (23) and (27) become 


f(t) = a Vt)t , 


(43) 


F(t) = H(t){— -(tot - — ) + e ^ 2 [(1 - 2 ^r-)sinQt + cosftt]} 

ait to to la 

U a) 

+ vVto ^[1 - ( sinfit + cosfit)e Vt ^ 2 ] . 


(44) 


Substitution of Eqs. (33) and (43) - (44) into Eqs. (40) and (41) 
yields : 


«/,n _ ■ r v/2 , H(x) r , N ^2, v V ■, vt/2 

A W - + f — " ~7T 11 " iT ‘ 2(v/2w) 1 - — > e c °sfir 


+ [<J> 0 + ~^2 H ( T ) " ^y]e VT ^ 2 s infix 
oto to 

+ ^[1 - 2(v/2oo) 2 ] + ^r- 

2o*a 


2uj 


(45) 


and 


r,/,s _ r v /2 . H(x) ri N ^2, v V vx/2 . „ 

BW - - {— ^ - -^-[1 - — - 2(v/2to) ] - - Z - T ^)e sxnfJx 


, n , v . vV. vx/2 
+ l4 , o + — 2 H ( x ' TjO® cosfix 

060 <0 


2fi to 


-^2 H(x) + 

OoJ OJ 


(46) 


where 


x H(x) = a (tp - tp ) > 0 
o 


(47) 


by Eq. (42) . 
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Upon substitution for A(i|0 » B(i{0» and f(t), F(t) , in 

accordance with Eqs . (45) - (46) and Eqs . (43) - (44), respectively, 
one obtains from Eqs. (30) - (32) the parametric solutions: 


^ — = oQip + ofix(<jj,t) - fttH(t) , 


(48) 


mm /oQe 


= H(t) - (1 - — )H(T)e" v(t T)/2 cosfi(t-i) 

1/a 


n 


- { ^"[(l + — )H(t) - 2aV] + (^) 2 afii|/ )e ^ x ^ 2 sinft(t-x) , 
ZU n sc o 

o 

(49) 


= H(t)[H(t) - e V(t T)/2 {(1 - 5-)H(T)cosn(t-T) 


n 


+ [(v/2fl)(l + — )H(x) - (v/ft)aV + (w/a) 2 crCty Isinn(t-x)} 
n o 

o 


. »_ Mii e-'' t/2 si„!!t)- 1 


n 

o 


(50) 


where : 


afix(iji,t) = fttH(t) - oQip - (vfi/co 2 ) [H(t) - aV] 

- {[1 - — - 2(v/2u) 2 ]H(t) + 2(v/2w) 2 aV - (v/2 0)afli|i }e _V (t_x) /2 sinft(t-x) 
n o 

o 

+ { (vfi/w) 2 [H(x) - aV] + }e T ^ 2 cos^(t-x) (51) 

o 

by Eq. (24). In Eqs. (48) - (51), x = (ip) is given by Eq. (47). 

Eq. (51) gives ip = tHx,t) for every point (x,t) of the wave. The 
boundary condition ip(x=0,t) in Eq. (37) is satisfied by Eq. (51) 
since t = x for x = 0 . 

Let nondimens ional variables and fields be introduced by the 


substitutions : 
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x/x -»■ X , t/t -*■ t , x/t ■+ x , \p/x \p , 
o o o o 


E/E -> E , v/V v , V/V + V , n/N -*■ n , 
o o o o 


(52) 


where 


x = l/oft , t = 1/Q , E = mw /aQe , V = 1/a , N - N . 
o o o o o 

(53) 

Further, if surface charges are absent (ty = 0) , the Eqs. (49) - 


become, for t > 0 and x > 0 : 


E(t,t) = ip + x(i)j,t) - t 
v(<p,t) = 1 - (l-a)e~ e(1 " e ) '^cos(t-*) 


(54) 


-6(1-B 2 ) h (l + a - 2V)e 


-eU-eV^t-ilO 


sin(t-ijj) , 


(55) 


n(^,t) = {1 - e 


-6(1-6 W -ip) 


[ (l-a)cos ( t— 1 ^») 


2 -Jj -1 

+ 6(l-6 2 ) 2 (1 + a - 2V)sin(t-iJ>) ] - a6(ij;)e ^ ^ t sint} 


where 


2 h 

x(<lst) = ( t-d>) - 26(1-6 ) (1-V) 


- [1 - a - 26 2 (1-V)]e 0(1 6 ) (t ^sin(t-ij>) 


(56) 


+ 26(1-6 ) 2 (1-V) e 


-6(1-6 ) " 2 ( t— xp) 


cos(t-i|>) , 


(57) 


and 


*l> = x , ip (x , t) >0 (x > 0) , 


(58) 


(51) 
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by Eq. (47). The dimensionless constants are defined as 

a = N/n , 6 = v/2u (59) 

o 

Eq. (57) gives ip = ijj(x,t) , where ifj(x,t) >0 by Eq. (58), for 

every point x of that region x > 0 which is occupied by electrons 

at time t . The limit, t^(x,t) = 0 , defines a function x = x(t) , 
which represents the moving position coordinate of the front of the 
electron gas. By Eqs . (54) - (56), the fields E(^,t) , v(ip,t) , and 

n(ip,t) depend exclusively on (t-i|/) , i.e. on x by Eq. (57).' The 

spatially periodic field configuration, which grows with the speed 
dx(t)/dt into the space x 2 0 , represents a (nonlinear) 
"neutralization wave". 

The Eqs. (54) - (58) are based on a nonlinear transformation which 
exists if J = 3x(ijj,t) /94/ ^ 0 or n(t/>,t) < 00 in Eq. (56), i.e. if 

0 < a < 2 + e , or: — < n /N < °° , 

2+e o 

2 _j. (60) 

e = exp[8(l-B ) 2 tt] - 1 « 1 , 

since the extrema of n(ifi,t) ^ occur when sin(t-r) = 0 . The 

position x(t) of the neutralization front is obtained from Eq. (57) 
as the limit x(t) = x(f -*■ 0, t) , 

x(t) = t - 2B(1-6 2 ) 3 s (1-V) 

- e~ B(1_e \ 2t {[l - a - 26 2 (1-V)]sint 

- 2B(l-6 2 )^(l-V)cost} > 0, t > 0 , 


( 61 ) 
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by Eq. (60). It follows for the speed of the neutralization front. 


= 1 - e 6(1 3 ) t [(l-a)cost + 6 (l-g 2 ) ^(1 + a - 2V)sint] 

dt ■ 


= v ^’ t) ^=0 


(62) 


In the region ahead of the neutralization front, x > x(t) , it 

is t < 0 and tb = 0 (ip =0) or i p = constant (ip $ 0) . 

o o 

Accordingly , 


n(x,t) = 0 , v(x,t) = 0 , x > x(t) , 

and 


(63) 


E(x,t) = E[x(t),t] + x - x(t) = x - t , x > x ( t) , (64) 


by Eq. (54), is the space charge field in the unneutralized ion space. 

In Figs. 2, 3, and 4, the nondimensional fields n(x,t) , E(x,t) , 
and v(x,t) [Eqs. (54) - (56)] of the neutralization wave are shown 

versus 0 < x £ x(t) at time t = lOn , with a = 0.1, 1.1 and 

2 1/3 -3 

V = 10 a. as parameters, and 6 = (m/M) / 4 tt = 1.107 x 10 for 

~ 2 1" 

mercury ions. Accordingly, it is x(t = 10 tt) = 10tt - 26(1-6 ) 2 * 

2 k 

(1-V) (1-exp [-6 (1-6 ) 2 10 tt]} = 10it . For other times t = t , the 

A A A 

fields of the neutralization wave exist only up to the point x = x(t) , 
e.g., the broken lines in Figs. 2-4 represent the neutralization front 
at time t = 2ir . The position coordinate x(t) of the neutralization 

front is given in dependence of t with a = 0.1, 1.1, and V/a = 

2 -3 

10 as parameters (6 = 1.107 x 10 ) in Fig. 5. 
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The momentum transfer between the electron gas and the ion beam 
[g-terms in Eqs. (54) - (56)] results in i) an amplitude asymmetry 
which is clearly visible in the case ct = 1.1 of Figs. 2-4, and ii) 
damping of the standing wave amplitudes as shown in Figs. 2-4. 

According to Eqs. (54) - (57), the wave fields exhibit for large x 
the symptotic behavior: 

n (x , t ) -* 1, v(x,t) -> 1 , 

2 1/2 

E(x, t) + -20(1-3 ) (1-V) , for x » x , (65) 

c 

where 

x c = (i-g 2 ) 1/2 /e . (66) 

is the nondimensional length of relaxation. Accordingly, the inter- 
component momentum exchange leads to a dissipation of the standing wave 
structure a few relaxation lengths x^ downstream of the injection 
plane (x = 0). This result means that a complete neutralization is 
achieved after a length Ax = 3x c for injection ratios a ^ = n Q /N ^ 1. 

The dimensional length of relaxation is by Eq. (53) 

x c = (l-g 2 ) 1/2 x o /6 = (2/ v) | i Q | /Ne . (67) 

For neutralization experiments, it is noted that the relaxation length 

x c decreases with increasing ion density N and decreasing electron 

current i Q . Note that the neutralization remains incomplete in the 

region 0 < x < x at all times if a = N/n ^ 1. 

— ~ c o 

In Figs. 6 and 7, the spatial attenuation of the nondimensional 

( 

amplitude of the electron density field n(x,t) is shown at times 
t ^ IO^tt for a = 0.1, 1.1, V = 10 2 a, and 3 = 1.107 x 10 [In the 
special case, a = N/n Q = 1> practically a relaxation free neutraliza- 
tion results within the wave, since n(x,t) = 1 by Eq . (56) for 
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0 < x < x(t) if gV << 1.] It is seen that the ion gas is essen- 

3 

tially neutralized at the point x = 10 it where n = 1. 

The reference values x , t , E , V , and N defined in Eq. (53) 

o o o o o 

determine the characteristic scales and actual magnitude of the fields 
of the neutralization wave. As a numerical example, consider a 
mercury ion beam and an electron injection current for which: 

N = 10^ cm 3 , i = 10 ^ amps cm 2 , g = 1.107 x 10 3 

Hence: 

-1/2 -3/2 

x = 3.691 x 10 4 (l-g 2 ) |i |N = 3.502 x 10~ 3 cm 
o 1 o 

-5 2 ~ 1 ^ 2 -1 / 2 -io 

t « 1.773 x 10 (1-6 ) N = 5.607 x 10 sec 

o 

1/2 - 1/2 

E q = 2.227 x 10 _4 (l-6 2 ) |i |n - 2.113 x 10 _2 cgsu = 6.339 x 10° volt cm' 1 . 


V = 2.082 x 10^ I i In 1 = 6.246 x 10^ cm sec 1 
o o 


9 -3 

N = N = 10 cm 
o 


and 


-3/2 

x = 3.691 x 10 4 6 -1 |i IN = 3.163 x 10° cm 
c 'o' 
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and v = 10 a showing spatial attenuation. 
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N/n = 0.1, 1.1 and V = 10 a showing spatial attenuation. 
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FIG. 5: Position coordinate x(t) of the neutralization front versus t for 
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FIG. 6: Spatial attenuation curve showing amplitude n(x^,t) of electron density wave at x^= 2irk, 

o 2 

k = 1, 2, 3, ... , and t ~ 10 tt for a = N/n = 0.1 and V = 10 a (neutralization: n 1). 





FIG. 7: Spatial attenuation curve showing amplitude n(x, ,t) of electron density wave at x = -r(2k+l) 
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V. INTERCOMPONENT MOMENTUM TRANSPORT AND ELECTRICAL 
CONDUCTIVITY OF COLLISIONLESS PLASMA 


ABSTRACT 

Based on the Lenard-Balescu equation, the interaction integral 
for the intercomponent momentum transfer in a two-component, collisionless 
plasma is evaluated in closed form. The distribution functions of the 
electrons and ions are represented in form of nonisothermal, displaced 
Maxwellians corresponding to the 5-moment-approximation. As an 
application, the transport of electrical current in an electric field 
is discussed for infrasonic up to sonic electron-ion drift velocities. 
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In a fully ionized plasma, the interaction radius of an electron 

or ion is of the order of magnitude of the Debye radius, D . In many 

plasmas of technical interest, such as glow discharges and high 

temperature, fully ionized gases, the Debye radius is large compared to 

- 1/3 

the mean particle distance, D >> n , where n is the electron or 
ion density. In this case, a charged particle interacts simultaneously 
with many particles of the plasma, since the sphere of interaction 
contains a large number N of electrons and ions, 

N ~ nD 3 » 1 . 

A kinetic equation which takes into consideration the many-particle 

1 2 

interactions in collisionless plasmas is that of Lenard and Balescu . 

1-2 

The formal mathematical deductions of this kinetic equation apparently 
give no meaningful explanations of its. physical foundations and its 
applicability to nonequilibrium situations in general. In the following, 
an application of the LB-equation is made which has the purpose of 

i 

determining the linear momentum exchange between the electron and ion 

1 

components and the associated electrical conductivity of a collisionless 

plasma. Since the LB-integro-differential equation can not be integrated 

in closed form because of its mathematical complexity, a simple plasma 

model is chosen by means of which intercomponent nonequilibrium effects 

can be studied as known from similar investigations based on the Boltzmann 
3-4 

equation . The electron and ion components are assumed to be in an 
approximate thermal equilibrium at their respective temperature 
(T r ^ T g ) in their individual center of mass system (<v r > f <v s >). 

The drift velocity between these components can be subsonic or 
supersonic. 
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where m and T rg are the reduced mass and temperature. This 
idealized theory should give information on the dependence of the 
intercomponent momentum transfer on the thermal nonequilibrium, 

T - T f 0 , and on the intercomponent drift velocity, <g > = <v > - <v > 

IT S ITS TC S 

At proper supersonic drift velocities, the collective interactions will 
become turbulent in the plasma (two-stream instability) . In this drift 
velocity region, the LB-equation is expected to break down or to give 
at best an approximate result with logarithmic accuracy. 
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THEORETIGAL PRINCIPLES 


The many particle collisions in a plasma can be regarded essentially 
as binary collisions complicated by the effects of the surrounding many- 
particle medium which becomes polarized by the respective interacting 
particle. If the time-dependent polarization of the plasma produced by 
this "test" particle is Fourier-analyzed, it is seen that the many 
particle collisions are equivalent to an interaction of the test 
particle with electrostatic waves (non-relativistic) . Evidently, these 
interactions occur at wave numbers k or impact parameters p (p k 
for which"* 


kD » 1 , p « D , 

where D is the Debye radius. The change of the distribution function 

f = f (v ) of particles of type "r" as a result of many-body 
it it r 

interactions with particles of type "s" is of the form 


3f 


/ s 


3t 


9v 


r/s 


( 1 ) 


j ^ is the flux of r-particles in velocity space which can be derived 
intuitively from the concept of a diffusion in velocity space**, 

+ oo + 00 f 9f f 3f . + + 

j , = (2e 2 e 2 /m ) j f f U- (— — - — ~^)dv , j . =0 (2) 

^r/s r s r m g £ « r £ s’ J s/s 

s r 


where 


^max kk6(ic*v - k*v ) 

r = jn . s d g 


k . k^ I c(k,aj) 1 2 
mm 1 ' ’ 


w=ic*v 


( 3 ) 
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e(k,to) is the longitudinal dielectric constant ( v^ << c) of the 


plasma 5, 


e£,«) = 1 + k' 2 l ^ /// (k-^)dv s . 

r,s 6 -tr ' oi-k-v^ 3Vg 


(4) 


Equation (1), together with Eqs. (2)- (4), represents the LB-equation 

for a fully ionized plasma consisting of an electron (r) and ion (s) 

component. [The sign tr in Eq. (4) indicates that in the integration 

the pole io-k*v =0 in the complex v -plane is to be traversed from 
p p 

below. ] 

The order of magnitude of the dielectric constant is estimated by 
evaluating the integral in Eq. (4) formally for a Maxwellian distribution 
of 3 -particles which gives 

k*v. 


1-2 


:(k,u> = k*v ) = 1 + l 


e 


(5) 


. 2 _ 2 ■f- ■+• 

r,s k D g k-v r -k*v e 


-2 2 2 

where D = 4-rm e /KT . Since the average expression <g> is of the order 
p p p p 

one for the equilibrium distribution and of the same magnitude for 

nonequilibrium distributions of g-particles, the dielectric constant 

can be replaced by unity except in the case where f is unstable"*, 

p 


e(k,ai = k*v ) = 1 , k^D \ i k^D^ >> 1 

r p 


( 6 ) 


The integration of Eq. (13) produces the (many-particle) interaction 
logarithm 


max ,, k 

dk , . max 


A " J ^=ln(^) = ln(DKT /|ee !) , 


(7) 


min 


min 


since 


k = D -1 , k = KT / | e e 
min max rs 1 r s 


( 8 ) 
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where 

D = (4tt l n^/KT )~ H . 

r,s e e 8 


( 9 ) 


In Eq. (8), k corresponds to the maximum impact parameter 
P m = D , i.e., ^ m -^ n represents an extrapolation of the inequality 

in Eq. (6) which is permitted in view of the weak, logarithmic 
k-dependence of A , Eq. (7). 

ITS 

— i 

By combining Eqs. (3) and (6), an approximation to the tensor U 

is obtained which will be used exclusively further on, 
k 


max 
^ /// 


k . 
min 



<5(k' 


v -k*v )dk 
r s 


( 10 ) 


The distribution functions of the charged particle species r and s 
are assumed to be displaced, nonisothermal Maxwellians (5 - moment 
approximation)^ ^ , 


f 

r 


x , r.3/2 
(v ) = n (— ) 
r r tt 



(v -<v >) 
r r 


f 

s 


f 

r=s 


(ID 


In addition, let some frequently used abbreviations be defined in 
advance^ ^ , 


m =mm/(m+m) 
rs r s r s 


T 

rs 


T T 
/ r . s. 
m ( — + — ) 
rs m m 
r s 


9 


a = m /2KT 
r r r 


a 

s 


a 

r=s 


a 

rs 


= a a / (a +a ) 
r s r s 


( 12 ) 


In connection with the analytical developments, the so-called error 

function is encountered which is defined by 

! x 2 

$(x) = 2ir 2 j e 
0 


dt 


(13) 



where 
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$(x) = j ^ 2 

2 * 2 v<,“* 


xe 


2n f C~l) m x 2la+ l 



X > 


( 14 ) 


m*=o 


(25^r , x < 2 


( 15 ) 
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INTERCOMPONENT MOMENTUM TRANSPORT 


The force density exerted by the s -component on the r-component 
through particle-wave interactions is according to Eqs. (1) and (2) 

+<» „ +» . . f 3f f 3f 


F . « -2e 2 e 2 ffl v t— *///t •(— — 1 - — — I )dv ]dv 
r/s r s JJJ r ^ m g £ m r £ s r 


(16) 


3v -CO 

r 


The 5 -moment-approximation to the distribution functions in Eq. (11) 
give 

f 3f f 3f a a a ^ ^ a . ^ 

2n n ( — r — ) [— (v — <v >) (v -<v >)] 

m m .->■ r s 2 l m r r ' m s s ' J 

s 3v r 9v it r s 

s r 


r -*■ \ 2 ,+ .2, 

exp[-a r (v r -<v r >) - a s ( v s - <v s > ) ] • 


(17) 


Let thermal velocities c , c , the relative velocity g , and 

r s * 3 & rs 

the reduced thermal velocity c rg be introduced by: 


c = v - <v > , c = v -<v> : 

r r r s s s 


(18) 


'rs 


v - v 
r s 


’ 8 sr 8 rs * 


-»■ 

c 


rs 


a a 

rs -*■ , rs -> 

c + c 

a r a s 


c = c 
sr rs 


(19) 

( 20 ) 


Hence, 


-*■ ->■ , rs.+ -> . 

c = c + (g - <g >) , 

r rs a °rs °rs 
r 


c = c - (g - <g >) , 

s rs a & rs B rs ' 


( 21 ) 

( 22 ) 


and 


c - c = g -<g> 
r s rs 6 rs 


( 23 ) 
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Since the Jacobians are J, = |[3(v ,v )/3(c ,c )]| =1 and 

x r s it s 

= |[3(c ,c )/3(c ,g )]| = 1 , the velocity space elements in the 
z r s rs rs 

various variables equal, 


dv dv = dc dc = dc dg 
r s r s rs rs 


(24) 


By means of Eqs. (18)-(22), the bracket expressions in Eq. (17) are 
transformed as 


a a a a ^ a 

r s -»■ , r s N -»- , rs -*• _ . 

— c c = ( ) c + (g - <g >) 


m r m s 
r s 


'm m rs m rs 
r s rs 


rs 


(25) 


and 


-*■2 -*-2 . N -»-2 , 2 
ac +ac =(a + a )c + a (g - <g >) 

rr ss r srs rs °rs to rs 


(26) 




/y ~y 


The relations c = c (v ,v ) and g = g (v ,v ) are given in 
rs rs r* s' °rs °rs r’ s 6 

Eqs. (20) and (19), respectively. Hence, the divergence operator becomes 


rs 3* + 3 ■ 


3v s 3c 
r rs 


(27) 


3g 


rs 


The tensor U in Eq. (10) is evaluated by means of spherical coordinates 
(1,0, <f>) with the polar axis parallel to g (k° = k/k) , 

ITS 

k _ 

max tt 2tt ^ 
dk rt to to 


U = — i— J ^ j j k° k° 6 (cos0) sin0d0d<j> 


J rs k . oo 
min 


-*■ -*■ 

SR. 


= iv(A /g ) / [^%in^8 + (2 cos^0 - sin^O)— - r - S -] 6(cos0)sin0 d0 , 
r s r s ' 


i.e. 


■+ -V 
2 • 2 

tt /. / \ /tt rs e rs N 

U = tt (A / g ) (6 - ) 

rs e rs 2 


3 rs 


g 


rs 


where A is defined in Eq. (7), and 6 . . = 1 or 0 for i = i 
rs ij 

or i i j . 


( 28 ) 
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The integral over v^ in Eq. (16) reduces to one over 

c r = v^_ - <v r > , since the term proportional to <v r > vanishes by 

Gauss's law. Since v r and. v g are independent variables, the 

operator can be put under the dv; -integral. Thus, one 

r s 

obtains, under consideration of Eqs. (17) and (25)— (28) , 

± , 22. . 2.3/2 . t“ f 

F . =-4ttc e (a a /it ) n n A J ...I 

r/s rsrs r s rs J J 

(c i 3- + -2— : )-<AB)]dJ dj i 

v rs a “rs a rs 6 rs 

r s 8 c 8 g 

rs °rs 


exp [-(a +a )c^ - a (g -<g >)^] , 

* r s rs rs °rs °rs 

, a r a s. , C rs ( *r 3 ' a ts ) - , , a r8, ( 8 rs -<S rs » + 

‘ 5T )[ I 3 ®rs + 3 8 rs 

r s “rs g rs rs g rs 


Accordingly, 


? , =-4ire^e 2 (a a /ir^J^^n n A (I.+f +I„ +I 9 ) 
r/s rsrs r s rs 1 12 21 2 


where 


a +°° 

*1 


— /--Jc — 7 — • (AB) dc dg , 

a s L J rs & rs rs 

rs 


X , 5 /// XjiTi -J- -(AS)dS rs ]dj rs - s . 


12 ' a a rs " ‘ 

r g -oo _oo 3c 


+°° +00 

/// J rs C//J 

-°° -OO 3g 

rs 


(AB)dg ]dc = 0 , 

6 rs rs * 
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H 


rs 


Now, 


cl SL 

-(AB) = 2A[-(a +a )c *B + (-^ - -^)-M 
8c r s rs m r m g 

rs 


and 


— — • (AB) = 2A{ -a (g -<g >) *B 

-> rs °rs rs 


9g 


rs 


a rs ( 8 rs‘ < 8 rs >) a r a S (g rs‘ C rs ) 

+ [__Es r s rs (-£-_«) rs rs ] } 

m J m m -i 

rs g r s g 


'rs 


rs 


by Eqs . (30) and (31). Since A is an even function of c 


rs 


follows 


a +°° (g *<g >) (g *c ) 

1 - -2a -ES- / .../ Ac rS 8 rs rs 

1 r m J i rs 


rs 


3 

’rs 


.->• -*■ . 

(c *<g >) 

v rs B rs , ,-*■ ,-y 

— ] dc dg 

g rs °rs 

& rs 


and 


a a +°° 

1 = 2 — E*-E* J .../ Aj (a [ 

2 am 1 ' ■' °rs rs 
r rs -“ 


.-*■ .2 2 

^ 8 rs* <8 rs >) <8 rs > 


3 

y 

’rs 


g 


rs 


(g »<2 >) 

, VB rs 8 rs i /*■ 

+ 5 } dc dg 

3 rs e rs 


g 


rs 


Upon performing the dc -integrations, the Eqs. (39) and (40) 

IT S 


to 


(36) 


(37) 


(38) 


it 


(39) 


(40) 

reduce 
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* „ 3/2 a rs /a s - a rs ( *rB' < ®rs >>2 

*1 ' - — 7 ~ ,3/2 >>l C 
rs (a +a ) -“ 

r s 


x l 


-y -y 

g '<8 > 

fa rs -»■ 

— g 


rs 


<g > 

- B rs 


g 


rs g 


]dg 


rs 


rs 


rs 


and 


3/2 a^ /a +°° -a (g -<g >)^ 

x_ 0 tt rsr trr rs 6 rs ®rs -*• 

*2 - 2 — 7 — 3/2 > >1 6 «rs 

rs (a +a ) -°° 

r s 


x {a [ 
rs 


r>- -*■ .2 

<8 rs’ <8 rs > 


J rs 


-*■ 2 ->-■»■ 

<g > e *<e > 

s rs , , B rs B rs , 

T~ 1 + __ 3 )d8 rs 


rs 


g 


rs 


( 41 ) 


(42) 


Let a spherical coordinate system (g ,a,g) with the polar axis 

rs 

parallel to < g rs > bn introduced, in which 

g rs = g rs (sina cos0, sina sing, cosa) , 

-y 2 

dg = g dg sina da d8 , 
rs rs rs 

and the substitutions. 


(43) 


t = cosa , 



(44) 


These operations transform the Eqs. (41) and (42) to: 


5/2 a /a 

-* _ 9 it rs s 

1 




m , A v 3/2 v 12 1(T s rs 

rs (a +a ) 
r s 


and 


5/2 a /a 
t „ 4 £ rs r 

2 m 


rs (a +a ) 
r s 


3/2 *- Y rs (G 23 G 21 ) + G 12 ]<8 rs > 


(45) 


(46) 
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where 


mn 


2 2 

+1 °° -(x +y -2y xt) 

r r 1 rs 1 rs m , n , 

J J e x dx t dx 

t=- 1 x=0 


(47) 


By means of an obvious substitution, this integral is brought into the 
form 


+1 Y 2 (t 2 -1) 

G - J H (t) e rS t n dr , 

mn J m 


(48) 


where 




-y T 
'rs 


e" 5 (C+Y rs T) m d? 


(4 9) 


1 • 6* j 


m 


°° r 2 

m-n r 


H X = y ( )(y x) j e 

m ' L n v n 'rs ; 

n=0 -y t 

rs 


/ e' C i 1 dK 


Hence , 


; 10 y 2 ®* Y rs? 


rs 


(50) 


Y~ - 1 


12 


—V 

rs tt * A/ \ , 1 rs 

~ 3~2 t <l' r s )+ — e 


rs 


' rs 


(51) 


21 


2y 2 - 1 h 1 -y 2 

2~ 2 *<*«> + — " 


2y 


rs 


rs 


(52) 


3 - 3y 2 ' + 2y 4 % 

rs rs it 


23 


2y 


rs 


i + , 3 

2y 


“ A 

Y_ - 3 -y 


rs 


2 

rs 


rs 


(53) 


and 


a 1 
s 1 

-y 

a 

r 2 


5/2 a 


= 2 


TT 


rs 


m , , ^3/2 ’rs 2 

rs (a r +a g ) 


-3 TT 2 r , , 

Y T l*(Y rs ) 


1 2 " Y rs, „ -»■ 


Y ~r e 
rs % 


1 < g > 

& rs 


( 54 ) 
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Substitution of Eq. (54) into Eq. (32) gives the following expression 
for the force density exerted by the s-component on the r-component: 


-1 


F , = - M(a <g > )t n m (<v > - <v >) 

r/s rs 1 e rs 1 ■ rs r r r s 


(55) 


where 


% - y 2 

M(y ) = y 3 [$(y ) ~ r Y e rS ] , 

rs 4 ’rs 1 u rs ; % 'rs 1 

7T 


(56) 


r 1 

rs 


= ir( 


c 2KT h n 
8 , rs , rs 


irm 


rs 


m 


n Q 
s rs 


e e 2 

JL ( JL± . A 

0 ' 17T ) 


*rs 2 v KT 


rs 


rs 


(57) 

(58) 


This result indicates that the intercomponent friction force 
^ r / s = (action-reaction principle) is in general a transcendental 

function of the drift velocity <g > = <v > - <v > [Eqs. (55 )— (56)] . 

ITS TO S 

The relaxation time t is an expression analogous in structure 

3-4 

to that for binary Coulo mb collisions [Eq. (57)]. Q is the many- 

rs 

body Coulomb transport cross section [Eqs. (58) and (7)]. 

The transcendental function M(y ) in Eq. (56) has the limiting 

Tl s 

properties [Eqs. (14)— (15) ] : 


M(Y rs > 

S 1 • Y rs 

« 1 , 

(59) 

M <Y rs ) 

- 3 % -3 

= ~r IT Y 
4 rs 

, Y >> 1 . 

rs 

(60) 


Accordingly j 


-1 


/ = r n m (<v > - <v >) , y « 1 

r/s rs r r r s ’ 'rs ’ 


(61) 


and 


_i <v > - <v > 

r 1 / 3 t t . r s 

r/s -T rs r m r^. 3/2^ ", + + 

4a <v > - <v > 

rs 1 r s 


T » Y » 1 . 
1 3 rs 


( 62 ) 
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It is recognized that the intercomponent friction force i) increases 

proportional to the drift velocity in the infrasonic region, |f . | ^ g 

it / s rs 

by Eq. (61), and ii) decreases inversely with the square of the drift 

y _2 

velocity in the suprasonic region, |F . | 'v g by Eq. (62). It 
should be noted that Eq. (62) is not an exact extrapolation to supra- 
sonic drift velocities, since the LB-equation is no longer rigorously 
applicable for supersonic drift velocities due to the appearance of 

the two-stream instability (y >1). The Eqs. (55) and (62) can 

r s 

be used, however, as approximations in the supersonic and suprasonic 

region of drift velocities, respectively. As a justification, it is 

noted that Eq. (55) and Eq. (62) give not only a qualitative but also 

a quantitative description (within the experimental uncertainties) of 

9 

the so-called run-away-ef feet . 


J 
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' / 

ELECTRICAL CONDUCTIVITY 

In order to discuss the concept of electrical conductivity in a 
simple way, a quasi-homogeneous steady-state electron-ion plasma 
(r = e , s = i) is considered in a frame of reference in which an 
electric field E exists but no magnetic field (B = $) . The electrical 
current density, 

i ■ l n 6 e g <V = - ne< g ei > , (63) 

e,i p p 

n.=n = n , e . = -e =e , 
l e ’ i e 

is related to the argument Y rs defined in Eq. (44) by 

Y ei = j/j T (64) 

where 

j = |j| » J T = ne/aj. = ne(2KT e /m e ) 1 { 2 (65) 

In the steady-state of a quasi-homogeneous plasma the intercomponent 
friction and electric force densities are in balance, F .. - neE = $ . 
The corresponding Ohm’s law is by Eqs. (55) and (63)— (64) : 


j M(j/j x ) = a E , j < j T , 


( 66 ) 


where 


a = (ne /m )t . 

e ei 


(67) 


and 


_ 2KT % 

x . = -r( ) n Q . , m << m. , 

ei 3 irm ' x ei ’ e i ’ 


( 68 ) 


^ei = A ei ■ A el ‘ « m, , 

e 


( 69 ) 
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by Eqs. (57) -(58) and (7). According to Eq. (66), the relation between 

current density j" and electric field E is nonlinear if. j is not 

small compared to the thermal current density [Eq. (56)] . 

In the infrasonic drift limit, j << , a has the meaning 

of an electrical conductivity [Eq. (67)]. If the electron drift 

.l' ~ t 

velocity is small compared to the thermal speed a". 2 = (2KT /m ) 2 , 

61 6 6 

Eq. (66) reduces to the familiar linear Ohm's law: 1 ^ 

j" = aE , j « j T . (70) 

It should be noted that Eq. (66) is not applicable to suprasonic 

drift velocities, • In the latter case, the run-away-effect 

occurs, and an Ohm's law does no longer exist in steady state. 

The results obtained in this investigation by means of the LB- 

equation justify the use of the Boltzmann equation with a Rutherford 

cross section in plasma kinetics as long as the dielectric constant 

e(k,o)) can be approximated by one. It can be shown that the collective 

effects are only important if the ratio of the electron to the ion 

2 

temperature is very large, > 10 . A similar conclusion concerning 

the significance of the collective contributions has been obtained in 

connection with the thermal relaxation in anisotropic two-temperature 

11 2 
plasmas. It is noted that electron temperatures, > 10 T^, are ex- 
tremely difficult to realize in fully ionized plasmas with stable dis- 


tribution functions. 
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